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ABSTRACT 

An  effective  data  interface  between  a  Victoreen  PIP-400 
pulse-height  analyzer  and  a  Hewlett-Packard  9  810A  calculator 
was  designed,  built,  and  tested.   A  calculator  program  was 
written  which  enabled  a  research  group  studying  laser 
propagation  in  the  marine  boundary  layer  to  conduct  rapid, 
local  processing  of  scintillation  and  extinction  data. 
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I.   INTRODUCTION 

The  potential  utility  of  laser  technology  in  the  areas 
of  line-of-sight  communication,  target  tracking,  and  wea- 
ponry makes  an  understanding  of  the  propagation  of  laser 
beams  near  the  air-ocean  interface  imperative.   The  effects 
of  a  turbulent  atmosphere  on  laser  transmission  were  theor- 
etically modelled  by  Tatarski  [Ref .  1]  et.  al.  and,  subse- 
quently, extensive  field  research  in  the  terrestrial  boundary 
layer  has  substantiated  the  bulk  of  the  theory.   Research 
data  in  the  marine  boundary  layer  has  been  scarce  and 
inconclusive. 

Atmospheric  effects  on  laser  propagation  may  be  divided 
between  laser-induced  phenomena  and  those  caused  by  meteoro- 
logical characteristics  independent  of  the  laser's  presence. 
The  first  classification  applies  only  to  lasers  of  sufficient 
power  to  heat  the  atmosphere  (thermal  blooming,  upwind 
drift) .   The  second  is  a  universal  grouping  affecting  the 
whole  electromagnetic  spectrum  regardless  of  power.   The 
atmospheric  effects  include  beam  refraction,  beam  spread 
(or  degradation  of  resolution) ,  absorption,  and  scattering. 
Refraction  may  be  further  categorized  according  to  the 
dimensions  of  the  atmospheric  whirls.   Large-scale  vortices 
cause  path  deviation  (beam  wander)  while  small-scale  vortices 
cause  time-intensity  variations  (scintillation)  and  degrada- 
tion of  resolution. 


The  Naval  Postgraduate  School,  Monterey,  has  been  involved 
over  the  last  four  years  in  a  continuing  effort  to  establish 
the  data  base  necessary  to  confirm  or  modify  existing  theory 
as  applied  to  the  air-ocean  interface.   Two  research  groups 
have  operated  simultaneously  but  independently;  the  optical 
research  group  and  the  atmospheric  turbulence  group  have 
attempted  to  correlate  determinations  of  CN  (the  refractive 
index  structure  constant)  through  the  respective  analysis 
of  light  intensity  distributions  and  meteorological  data. 
The  optical  research  group  has  concurrently  investigated 
scintillation,  extinction  (scattering  and  absorption)  and 
the  atmospheric  modulation  transfer  function  (a  measure 
of  beam  spread  and  wander) .   Theses  on  the  effects  of  aperture 
averaging,  determination  of  an  atmospheric  extinction  con- 
stant, and  determination  of  an  atmospheric  MTF  for  a  black- 
body  source  were  completed  over  the  same  period  as  the 
author's  [Refs.  10,  11,  12]. 

A  major  problem  for  the  optical  group  was  the  inordinate 
amount  of  time  necessary  to  process  scintillation  and  extinc- 
tion data.   The  methods  employed  depended  heavily  on  operator 
analysis  of  raw  data  or  extensive  manual  input  to  computers. 
While  the  processing  theory  being  applied  was  appropriate, 
the  time  penalty  and  the  inherent  high  probability  of  error 
had  to  be  minimized.   The  author  concentrated  his  efforts 
in  this  area  employing  only  the  basic  results  of  laser 
propagation  theory. 


II.   PROJECT  BACKGROUND 

Various  members  of  the  Naval  Postgraduate  School's 
faculty  have  been  conducting  research  on  atmospheric  pro- 
pagation phenomena  since  1972.   The  first  experiments  inves- 
tigated the  atmospheric  MTF  and  scintillation  in  the  basement 
corridor  of  Spanagel  Hall.   In  the  scintillation  experiment 
a  laser  beam  reflected  from  a  flat  mirror  placed  at  the 
opposite  end  of  the  137  meter  passageway  and  entered  a  pin 
diode  detector  located  near  the  laser  source.   The  detector 
output,  proportional  to  the  incident  beam  intensity,  was 
passed  through  a  log- amplifier  and  recorded  on  magnetic  tape. 
The  data  was  later  periodically  sampled,  digitized,  and 
accumulated  by  a  pulse-height  analyzer;  the  resultant  dis- 
tribution was  Gaussian,  as  theory  had  predicted.   Refractive 
index  structure  constants  were  calculated  from  photographs 
of  the  distributions  by  measuring  the  curves'  full  width  at 
half  maximum  height.   The  MTF  measurement  used  a  six  inch 
parabolic  mirror  mounted  as  a  Newtonian  telescope.   The 
image  from  a  laser  at  the  other  end  of  the  passageway  was 
scanned  across  a  slit  using  a  cam-driven  rocking  diagonal 
mirror.   The  output  of  a  photocell  was  fed  to  a  boxcar  inte- 
grator to  generate  time-averaged  data,  which  was  then  recorded, 
Turbulence  in  the  hall  was  somewhat  controllable  and  detec- 
table reduction  of  MTF  with  turbulence  was  observed. 
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In  late  1972  a  group  of  interested  researchers  discussed 
possible  environmental  experiments.   The  first  such  experi- 
ment was  completed  10  January  1973  by  faculty ,  staff,  and 
students  and  involved  shining  a  low-power  (mW)  red  laser 

o 

(6328  A)  across  a  4.3  kilometer  range  at  the  southern  end 
of  Monterey  Bay.   The  research  vessel  Acania  was  anchored 
at  the  approximate  midpoint  of  the  beam  path  to  take  atmos- 
pheric turbulence  measurements  with  onboard  sensors.   In  a 
second  experiment  at  San  Nicholas  Island  (3  March  1973)  low 
power  laser  beams  were  transmitted  from  the  northwest  tip 
of  the  island  to  the  Acania  anchored  offshore.   The  qualified 
success  of  these  experiments  completed  with  equipment  already 
owned  by  the  school  was  encouraging  enough  to  prompt  the 
submission  of  a  research  proposal  to  the  Naval  Ordnance 
Laboratory  (White  Oak,  Md.)  for  fiscal  year  1974. 

Experimental  apparatus  has  changed  significantly  since 
the  Navy  lab  started  sponsoring  the  research.   A  Cassegrainian 
telescope  system  with  an  eighteen  inch  parabolic  mirror  and 
an  elastically  mounted  scanning  mirror  (driven  electromag- 
netically  by  a  shaped  waveform)  was  used  in  the  most  recent 
MTF  measurements.   An  enclosed  trailer  was  built  to  transport 
and  protect  the  telescope  and  associated  equipment;  aluminum 
support  legs  that  extended  to  rest  directly  on  the  ground 
negated  the  previously  persistent  problem  of  telescope  motion 
due  to  operator  movement.   More  recently  a  bus  has  been  used 
which  holds  most  of  the  on-shore  optical  and  meteorological 


11 


equipment  and  employs  the  same  independent  telescope  suspen- 
sion idea.   Ship  motion  originally  made  it  difficult  to  keep 
even  a  diverged  laser  beam  on  the  telescope  long  enough  to 
collect  maningful  MTF  data.   A  gyro-stabilized  tracking  system 
was  designed  to  go  aboard  the  Acania.   At  first  it  was  manned 
by  two  operators  who  sight-aimed  a  gyro-mounted  laser  with 
hand-operated  servo  joy  sticks;  one  controlled  horizontal 
and  the  other  the  vertical  train.   Little  useable  data  was 
obtained  with  this  system  because  of  problems  with  human 
coordination  and  fatigue.   The  solution  was  a  quadrant  detec- 
tor servo  system  which  will  lock  on  a  laser  beam  once  the 
operator  has  aligned  it  with  the  target.   Shipboard  scin- 
tillation detectors  are  also  mounted  on  the  gyro  frame  and 
thereby  maintain  a  constant  relative  position  in  the  incident 
beam  path.   The  gyro  has  been  housed  in  an  aluminum  shed, 
when  aboard  the  Acania,  to  protect  both  equipment  and  operator 
from  the  elements . 

Scintillation  data  was  best  taken  aboard  the  Acania;  a 
detector  with  a  wide  angle  of  acceptance  could  continuously 
sample  a  diverged  beam  from  ashore  more  readily  than  the 
reverse  because  of  the  effect  ship  motion  would  have  on  the 
direction  of  propagation  of  a  beam  originating  on  board. 
The  single-detector  system  originally  employed  required 
excessive  attention  during  data  runs  because  light  from 
environmental  sources  refracts  into  the  beam  path  and  can't 
be  distinguished  from  that  of  the  monitored  laser.   Since  all 
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The  bus  with 
telescope  in 
position.  Transit 
to  the  right  was 
used  in  range-to- 
ship  calculation. 
In  background  stands 
meteorological  tower 
supporting  sensors. 


Naval  Postgraduate 
Research  Vessel  Acania 
with  meteorology  gear 
deployed.   The 
aluminum  shed  which 
shelters  the  gyro 
is  aft  of  the 
pilot  house. 
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Right  side  of 
telescope  showing 
fan-scan  unit, 
laser,  detectors, 
tracking  unit, 
and  associated 
electronic  gear. 


Gyro-stabilized 
tracking  unit.   The 
detectors  are 
mounted  on  top  and 
the  lasers  are 
hidden  from  view 
below. 
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Oscilloscopes , 
Nuclear  Data  Analyzer, 
log-amplifier,  and 
demodulator;  here 
shown  set  up  for 
meteorological  data 
processing  on  the 
bus . 


Dry  lab  on  the 
Acania  with 
scintillation/ 
extinction  data 
processing  gear: 
PIP-400,  PAR 
amplifiers, 
demodulator,  log- 
converter,  HP9810A 
calculator  and 
plotter,  Oscillo- 
scope &  freq 
analysis  tape 
recorder . 
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photocell  data  passed  through  a  log- amplifier  a  skewed 
distribution  would  result.   In  an  attempt  to  negate  this 
effect  the  beam  had  to  be  interrupted  ten  seconds  out  of 
every  minute  to  zero-out  detector  voltage  generated  by  back- 
ground radiation.   A  dual-detector  system  was  built  with  the 
idea  that  the  output  of  one  pin  diode  detecting  only  background 
could  be  subtracted  continuously  from  that  of  another  unit 
aimed  at  the  laser  beam. 

Since  the  detectors  (supposedly  identical)  produced 
different  proportional  outputs  extensive  voltage  biasing 
was  necessary.   An  even  greater  problem  was  that  ship  motion 
continuously  changed  the  position  of  the  detector  relative 
to  the  beam  center.   Subsequently  another  single-detector 
system  was  designed,  this  one  employing  a  chopped  laser  source 
and  a  demodulator.   The  chopping  was  accomplished  at  the  source 
by  reflecting  the  laser  light  first  off  of  a  cylindrical 
mirror  which  fanned  the  beam  horizontally  and  then  off  of 
a  flat  mirror  vibrating  at  approximately  two  kilohertz  which 
made  the  beam  scan  vertically.   The  demodulator  unit  took 
the  detector  output  and  electronically  subtracted  the  minima 
(background)  from  the  maxima  (laser  &  background)  and  sent 
the  difference  to  the  log-converter.   Detector  sensitivity 
was  recently  increased  by  replacing  the  pin  diode  with  an 
avalanche  detector.   An  FET  follower  converts  the  high- 
impedance  output  of  the  detector  to  a  low  impedance  which 
can  be  transmitted  a  greater  distance  through  coaxial  cable. 
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This  has  allowed  all  data  processing  gear  on  the  Acania  to 
be  positioned  below-decks  in  the  dry  lab. 

Processing  of  scintillation  data  was  initially  best 
accomplished  through  physical  measurement  of  photographs 
of  resultant  distributions  displayed  on  the  PIP-400  pulse 
height  analyzer.   Attempts  to  record  and  subsequently  analyze 
data  with  integral  methods  on  a  computer  suffered  from  poor 
digitizing.   Digital  output  could  be  obtained  directly  from 
the  PIP-400  by  a  fast  printer  or  from  a  Nuclear  Data  128- 
channel  analyzer  by  teletype.   An  attempt  was  made  to  fit  a 
Gaussian  curve  to  the  data  using  least-squares  method.   This 
analysis  proved  highly  successful  when  extraneous  voltages 
did  not  skew  the  experimental  distribution  but  was  painfully 
slow  because  data  had  to  be  manually  input  to  the  computer 
point  by  point.   The  most  recent  improvement  has  been  the 
direct  interfacing  of  the  PIP-400  with  a  compact  computer 
programmed  to  process  scintillation  and  extinction  data  on- 
scene  within  six  minutes  of  sampling.   This  was  completed 
by  the  author  and  is  the  subject  of  this  thesis. 
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III.   THEORY 

A.   SCINTILLATION 

The  theory  of  optical  propagation  in  a  turbulent  atmos- 
phere has  been  extensively  documented  and  those  portions 
employed  by  the  Naval  Postgraduate  School  optical  propaga- 
tion group  have  not  changed  significantly  since  the  project 
began.   A  thesis  submitted  by  B.  C.  Haagensen  [Ref.  3]  in 
June,  1973  contained  an  excellent  outline  of  the  pertinent 
theoretical  derivations  and  an  extensive  bibliography.   The 
author  does  not  intend  to  list  or  verify  familiar  theoretical 
developments  but  rather  to  present  the  significant  end 
products  vital  to  data  processing. 

Theory  states  that  atmospheric  turbulence  is  a  random, 
non-linear,  three-dimensional  phenomenon.   A  turbule  is 
ideally  assumed  to  be  an  isotropic,  homogeneous  pocket  or 
swirl  of  air  with  a  particular  index  of  refraction.   The 
inertial  subrange  delineates  a  region  of  turbule  size  in 
which  refraction  is  theoretically  described;  it  is  bounded 
on  one  side  by  turbules  so  small  (a  few  millimeters)  that 
viscous  dissipation  of  energy  causes  their  immediate  collapse 
and  on  the  other  by  eddies  so  large  (order  of  meters)  that 
they  are  no  longer  homogeneous  and  energy  is  added  to  the 
turbulence  spectrum.   In  the  inertial  subrange  large  turbules 
form  and  transfer  energy  to  successively  smaller  turbules 
with  minimal  losses  through  viscous  dissipation. 
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Tatarski  [Ref.  1]  described  the  three-dimensional  power 
spectral  density  of  turbulence,  $  (K) ,  with  the  equation: 

2 
*n(K)  =  0.033  CN2  K"11/3  exp  [^] 

K 

m 
where : 


CL,  =  refractive  index  structure  constant 
K  =  spatial  wave  number  =  — 

I     =  coherence  length  of  a  turbule 

2-tt 
K   =  -r— ;  l~   =  lower  limit  of  inertial  subrange 
m   lQ        ^0  y 


CN  obviously  parameterizes  the  total  amount  of  energy 
in  the  turbulence.   The  random  distribution  of  turbulence 
cannot  be  forecast  and  must  be  described  statistically.   A 
beam  is  modulated  in  a  multiplicative  manner  when  it  passes 
through  successive  turbules  of  differing  index  of  refraction 
It  follows  that  the  natural  log  of  the  beam  intensity  is 
altered  additively  and  should  follow  a  normal  probability 
distribution  as  a  result  of  the  central  limit  theorem.   C., 
is  directly  related  to  the  structure  function  of  the  refrac- 
tive index  (the  mean  square  difference  of  the  index  of 
refraction  at  two  points  a  distance  d  apart)  by: 


Dn(d)  =  <[n(r)  -  n(r+d)]2>  =  CN2  d2/3 
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After  lengthy  derivation  and  a  number  of  assumptions  [Ref. 
4,  pg.  10]/  CL,  can  be  seen  to  be  directly  dependent  on  the 
variance  of  the  log  intensity  distribution  by: 


where 


CN  =  1.42  K"7'/12  a"11'12  C 


2  2 

Co    =  <[ln(I)  -  ln<I>]  >  =  variance 


=  optical  propagation  path  length 


CLj  can  also  be  calculated  by  a  statistical  determination 

2 

of  the  temperature  structure  constant,  C   ,  since  tempera- 
ture is  the  dominant  factor  in  a  meteorological  formulation 
of  the  index  of  refraction: 


where 


CN  =  (79.0  ^  x  10  6)  Cr 


P   =  mean  pressure  in  millibars 


T   =  mean  temperature  in  K 


2    <[T(r)  -  T(r+d) ]2> 
CT   =       d2/3 
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A  major  objective  of  the  project  is  to  study  the 
correlation  between  C  values  obtained  optically  and 
meteorologically . 

B.   GAUSSIAN  MODELLING 

According  to  theory,  periodic  sampling  of  a  laser  beam 
propagating  in  the  inertial  subrange  should  approximate  a 
lognormal  distribution;  if  the  natural  log  of  the  detector 
voltage  (proportional  to  beam  intensity)  defines  abscissa 
values  and  the  ordinate  marks  the  number  of  samples  returning 
a  specific  voltage  the  resultant  plot  should  be  of  the  form: 


t  ~\2 

y  =  C  exp  [-1*  ~  x)    ] 

2a 


where : 

C   is  a  constant 

x   is  natural  log  of  detector  voltage 

x  is  the  average  natural  log  of  detector  voltage 

2  2 

a      is  the  variance  (C«T  ) 

Taking  the  natural  log  of  each  side: 

2  -2       7 

In  y  =  -  iL^.  +  *|  +  (inc 2j)    =   ax  +  bx  +  c 

2a  a  2a 
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where : 


a  = j 

2a 


b  -  — 


-2 
c  =  In  C  -  — 


2a2 


Experimentally  the  natural  log  of  detector  voltage  axis 
is  divided  by  the  PIP-4  00  into  two  hundred  channels  of  equal 
width  and  the  y  axis  is  the  number  of  counts  (or  times)  the 
equivalent  detector  voltage  (or  beam  intensity)  was  detected 
over  the  sampling  period.  The  number  of  counts  in  each  chan- 
nel after  a  data  run  can  be  assumed  to  be  Poisson-distributed; 
for  one  measurement  the  68%  confidence  interval  is  y  ±  /y. 
Since  the  natural  log  of  y  is  taken  to  obtain  the  quadratic 
form,  the  uncertainty  in  ln(y)  becomes: 


d[ln(y)]  =  Si  =  G.   =  ± 

y    y    /y 


So,  for  the  log  distribution,  the  channels  with  the  largest 
numbers  of  counts  have  the  narrowest  confidence  intervals 
(smallest  uncertainty)  and  should  be  weighted  more  heavily. 
The  effect  of  narrowing  confidence  intervals  with  increasing 
counts  is  easily  illustrated  by  plotting  them  on  log  paper. 
In  the  original  program  it  had  been  decided  to  weight  the  data 
inversely  with  the  uncertainty,  i.e.  multiply  by  /y. 
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The  ln(y)  distribution  resulting  experimentally  can  be 
approximated  with  a  polynomial  by  means  of  least-squares 
fitting.   The  method  seeks  to  minimize 


where: 
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E   w.[ln(y.)  -  ax.2  -  bx.  -  c]2 
i=0   x  1      x 


x.  =  channel  number  i 


y .  =  counts  in  channel  i 

1 l 


w.  =  weighting  of  data  in  channel  i,  /y~ 


Setting  the  partial  derivatives  with  respect  to  a,  b,  and 
c  equal  to  zero  and  cancelling  common  terms  gives  three 
equations  with  three  unknowns : 


2  4  3  2 

Ew.x.    lny .    =   aEw.x.      +  bEw.x.      +   cEw.x. 
l    l        -*  i  ii  ii  ii 


3  2 

Ew.x. lny.    =  aEw.x.      +  bEw.x.      +   cZw.x. 
iiJi  ii  ii  ii 


2 

Zw.lny.    =  aEw.x.      +   bEw.x.    +   cEw. 
iJi  ii  ii  i 


Applying  Cramer's  Rule  [Ref.  6]  to  the  system  determines 
the  coefficients: 
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det(A1) 
a  = 


det(A) 


det(A9) 
b  = 


det(A) 


det(A3) 
c  = 


det(A) 


where 


4     2 
det(A)   =         Ewixi  ^wj_xi  {Ewi-Ewixi> 


3  2  3 

+    Ew.x.     {Ew.x.    Ew.x. -Ew.x.    Ew.} 
i   l    L       iiiiiiij 


2  3  2    2 

+   Ew.x.     {Ew.x.    Ew.x. -(Ew.x.    )     } 
i    l    L       l    l         l    l  i    i         J 


3 

det(A, )    =   Ew.x.     {Ew. x. Ew. lny . -Ew. Ew. x. lny . } 
1  1111      i      -1 1        i      l    l      JiJ 


2  2  2 

+Ew.x.     {Ew. x. Ew. x. lny . -Ew. x.    Ew.x.    lny.} 
iiu      li      li      -'iii        ii        Jr 


2  2 

+    Ew.x.  {Ew.x.     Ew.  x.  lny . -Ew.  x.  Ew.  x.    lay  .  } 
iiliiii      -'i         i    i      i    i      HiiJ 


4 
det(A9)    =    Ew.x.     {Ew.  Ew.x.  lny  .-Ew.x.  Ew.  lny.  } 

to  X  X  XXX  X  J_  J_  -L  J_ 


3  2  2 

+   Ew.x.     {Ew.x.     Ew. lny .-Ew. Ew.x.     lny.} 
1111         iJi         ill         2 1 


2  2  2 

+   Ew.x.     {Ew.x. Ew.x.    lny. -Ew.x.    Ew.x. lny.} 
ii  ii      ii         2  i         ii         ii      J  l 
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4  2 

det(A-)    =   Ew.x.     {Ew.x.    Ew. lny . -Ew. x. Ew. x. 
3  1111         iJi         1111 


+Zwixi  {ZwixiEwixi2inyi-Zw.x. 3Ew.lny. } 


+  Ew.x.  {Ew.x. 3Ew.x. lny. -Ew.x. 2Ew.x.2lny. } 

11      11     11   ■'l     11     11    "*! 


Exactly  the  same  forms  for  the  coefficients  result  from 
applying  Gauss- Jordan  elimination  [Ref .  6] .   A  total  of 
eight  distinct  weighted-data  summations  form  the  basic  ele- 
ments which  combine  to  give  the  coefficients.   Coefficient 

"a"  and  a  channel  number-to- log  voltage  conversion  factor 

2 
deliver  the  lognormal  variance  (C^T  )  which,  with  the  optical 

wavelength  and  range,  determines  the  refractive  index  struc- 
ture constant  (CN) .   Coefficients  "a",  "b" ,  and  the  conversion 
factor  fix  the  average  log  voltage  value  (In 

C.   EXTINCTION 

Extinction  is  the  loss  of  time-averaged  intensity  with 
distance  by  a  beam  propagating  through  the  atmosphere.   No 
extinction  measurements  had  been  taken  by  the  optical  propa- 
gation group  before  September,  197  5.   Assuming  an  exponential 
dependence  and  spherical  spreading,  extinction  can  be 
represented  by  an  equation  of  the  form: 

xoRo2 
ir  u 

where: 
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I   =  intensity  at  distance  R 

Iq  =  intensity  at  distance  R- 

RQ  =  some  arbitrary  standard  distance 

R  =  some  arbitrary  distance 

y   =  the  extinction  coefficient 


An  equivalent  form  of  the  equation  is 


I  =  -^[I0R02exp(yR0) ]  exp(-yR) 
R 


If  R0  is  some  small  distance  (on  the  order  of  a  meter) ,  the 
exponential  with  RQ  approaches  one  and  the  above  form  can 
be  simplified  to: 


C 
I  =  — yexp(-yR) 

Rz 


where 


2 
C  =  IQRn      (a  measure  of  source  intensity) 


Taking  the  natural  log  of  both  sides: 


lnl  =  InC  -  21nR  -  yR 


Solving  for  y : 
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In  (-£*•) 

InC  -  21nR  -  lnl      IR 


R  R 

The  constant  C  must  be  determined  at  the  time  of  the  experi- 
ment and  the  time-averaged  intensity  and  range  must  be 
measured.   Since  detectors  generate  voltages  proportional  to 
the  incident  intensities,  voltage  can  be  substituted  for 
intensity  with  meaningful  results.   Phil  Parrish  investigated 
the  extinction  constant  (C)  and  coefficient  (u)  in  a  con- 
current thesis  [Ref .  10] . 

A  statistical  analysis  was  done  to  determine  the  extinc- 
tion coefficient's  degree  of  sensitivity  to  error  in  C, 
R,  and  I : 


M(c,R,r)  =  lnc  -  21f  -  lnI 


dy  "  JC   dC  +  3R  dR  +  31  dI 


IB     -    if  (dV\   2      An2  (^^\   2      AT,2     4.       ($JtL\  2      AT2     1    1/2 

T  -   y{ ( 3cJ       dC      +    ( 3R}       dR     +    ( Jl]       dI     } 


Checking  sensitivity  to  error  with  the  true  value  of  y 

-3 
assumed  to  be  1  x  10   : 
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Obviously,  errors  in  all  quantities  have  more  effect 
at  shorter  ranges.   Errors  in  the  extinction  coefficient 
appear  to  depend  linearly  on  independent  errors  in  the 
extinction  constant,  time-average  intensity,  and  range  but 
decrease  with  increasing  true  range.   Errors  in  range  have 
a  much  greater  effect  than  equivalent  errors  in  either  of 
the  other  quantities.   The  last  two  calculations  in  the 
previous  table  use  maximum  error  values  estimated  to  be 
appropriate  with  present  equipment  and  procedure. 


31 


IV.   INTERFACING 

The  fitting  of  normal  distributions  to  scintillation 
data  through  a  computer  program  using  least-squares  method 
had  proven  highly  successful  but  also  extremely  slow.   The 
PIP-400  pulse  height  analyzer  used  to  accumulate  data  dis- 
tributions had  the  design  capability  to  output  the  information 
through  a  50-pin  plug  in  parallel  BCD  (binary  coded  decimal) 
readout.   The  HP9  810A  calculator  in  which  the  processing  took 
place  had  an  optional  parallel  input  interface  card  which 
accepted  BCD  input  through  a  50-pin  plug  at  the  same  voltage 
levels  as  put  out  by  the  PIP-400.   It  was  thought  that  a 
simple  interface  matching  the  proper  lines  would  obviate  the 
pokey  manual  input  of  data.   A  number  of  hidden  problems 
were  found  which  made  the  project  much  more  difficult  than 
expected. 

The  first  problem  was  that  both  machines  had  been  designed 
as  masters.   The  analyzer  sent  a  high-to-low  (voltage)  print 
command  and  waited  for  a  low-to-high  flag  from  the  peripheral 
(which  told  it  to  step  to  the  next  channel) .   The  calculator 
sent  a  high-to-low  voltage  command  to  an  external  source 
asking  for  data  and  waited  for  a  low-to-high  flag  which  told 
it  the  information  was  ready.   The  solution  was  to  invert 
both  sets  of  signals  so  that  the  data  request  from  the  calcu- 
lator appeared  as  a  peripheral  flag  to  the  analyzer  and  a 
data  print  command  from  the  PIP-400  appeared  as  an  information 
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ready  flag  to  the  HP9810.   If  the  analyzer  function  knob 
is  set  to  "PRINT"  before  the  calculator  asks  for  data  the 
first  channel  will  be  ignored,  but  since  the  first  (and 
last)  five  channels  of  the  distribution  are  zeroed  by  the 
program  anyway  (because  of  the  high  probability  of  extrane- 
ous counts  in  these  channels)  the  deletion  has  negligible 
effect.   Since  the  voltage  signal  levels  for  both  machines 
matched,  the  modification  was  completed  with  a  single  hex 
inverter  chip  (Signetics  7404)  using  two  of  the  six  inverters 
available.   A  six  volt  supply  lead  from  the  analyzer  and  a 
100-ohm  resistor  filled  the  inverter's  power  requirement 
(5  ±  .25  volts) . 

Matching  the  data  lines  proved  more  of  a  problem.   Data 
from  the  analyzer  in  parallel  readout  consists  of  six  data 
digits  followed  by  three  address  digits.   Digital  information 
is  relayed  in  BCD  which  requires  four  bitlines  for  each 
decimal  digit.   The  PIP-400  puts  out  a  low  current,  high 
impedance  data  signal  while  the  HP9810  uses  a  pull-up  resis- 
tor to  bring  an  expected  high  current,  low  impedance  signal 
up  to  five  volts  for  a  "1"  bit.   The  calculator  reads  any- 
thing over  2.8  volts  as  a  "1"  bit  and  anything  less  as  a 
"0"  bit.   With  the  command  lines  being  inverted  by  the  chip 
and  the  thirty- four  information  bit  lines  properly  mated 
(inside  a  junction  box  with  two  50-pin  plugs)  a  calculator 
print-out  of  analyzer  data  showed  all  9's;  this  meant  that 
the  pull-up  resistor  in  the  calculator  was  raising  all 
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analyzer  output  signals  to  over  2.8  volts.   This  supposition 
was  substantiated  when  the  correct  data  was  properly  read 
out  of  the  analyzer  by  a  fast  printer  through  the  same  output 
plug. 

The  first  solution  attempted  was  to  install  a  grounded 
quarter-watt  resistor  in  parallel  with  each  signal  line  to 
pull  the  voltage  of  a  "0"  bit  below  2.8  volts.   The  address 
numbers  were  successfully  printed  using  a  5600  ohm  resistor 
on  each  lead  but  the  data  digits  were  decidedly  non-uniform 
in  their  individual  bit  resistor  requirements.   Comparison 
of  calculator  printouts  with  fast  printer  readouts  were  used 
to  determine  the  necessary  resistance;  lead  requirements 
varied  between  4  700  and  8200  ohms  for  the  calculator  to  print 
data  digits  accurately.   Variable  (0  -  10,000  ohm)  poten- 
tiometers were  used  to  find  the  range  of  resistance  values 
between  those  which  fixed  a  bit's  output  at  either  "0"  or 
"1" .   The  proper  resistance  span  was  smaller,  in  some  cases, 
than  that  separating  successive  standard  resistors  and  was 
found  to  vary  with  time  (possibly  temperature  dependent) . 
Upon  testing  line  voltage  levels  it  was  found  that  while  the 
resistors  brought  a  "0"  bit  from  3.0  to  2.7  volts  it  also 
brought  a  " 1"  bit  down  from  5.0  to  2.9  volts  (approximate 
values) .   Since  the  fixed-resistor  interface  was  obviously 
too  finicky  it  was  abandoned. 

The  next  scheme  attempted  employed  operational  amplifiers 
(Signetics  yA747)  constructed  on  silicon  chips  as  unity-gain 
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voltage  followers  and  succeeded  in  generating  the  high 
current,  low  impedance  input  the  calculator  was  designed 
to  use.   Seventeen  op-amps  were  required  to  handle  the  thirty- 
four  bits  of  information  which  made  up  each  data  point.   The 
chips  were  all  connected  in  parallel  with  a  100  mA  dual  power 
supply  (run  on  115  volt  AC)  which  generated  the  ±  15  volts 
direct  current  necessary. 


input  A        k 

^— ^— ■  ii  r  li 


input  B 


'1       14  • 


■•2      13  •" 


••3       12  •- 
♦4       11* 


♦5      10  •- 

-•6       9  •- 


♦7        8  * 


#>    output  A 


V  =  -15   Volts 


operational  anplif  ier 
socket  wiring  for  use 
as  a  voltage  follower 


(A  &  B  are  BCD  bit  lines) 
output  B 


V  =  +15  Volts 


With  the  op-amps  and  additional  wiring  the  circuitry  was 

too  bulky  to  be  contained  in  the  original  junction  box;  after 
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the  system  had  proven  itself  in  the  field  it  was  redesigned 
and  built  with  a  junction  box  large  enough  to  hold  all  the 
components  internally. 

The  interfacing  of  analyzer  and  calculator  cut  the  time 
to  input  each  data  point  from  seconds  to  milli-seconds; 
with  a  200-channel  analyzer  (PIP-400)  this  meant  a  savings 
of  five  to  ten  minutes  processing  time  on  each  run  and  made 
timely  on-scene  analysis  of  experimental  distributions 
possible. 
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The  remnants  of  the  first  interface  using  operational 
amplifiers.   While  successful  it  was  too  vulnerable  in 
an  operational  environment. 


The  interface  in  its  final  form.   To  the  right  are  the  op- 
amps  and  inverter  chip,  to  the  left  is  the  power  supply. 
50-pin  plug  on  right  goes  to  calculator,  one  on  left  to 
analyzer;  outlet  plug  is  for  dual  power  supply. 
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V.   COMPUTER  PROGRAM  DEVELOPMENT 

A.   MODIFICATION  OF  THE  ORIGINAL  PROGRAM 

While  mechanically  interfacing  the  pulse  height  analyzer 
(PIP-400)  and  the  computer  (HP9810)  saved  a  great  amount  of 
time  previously  used  to  input  data,  the  computer  program 
which  processed  the  data  was  itself  inefficient.   An  inor- 
dinate amount  of  time  was  spent  labelling,  giving  operator 
instructions,  and  searching  for  subroutines.   The  program 
had  been  written  in  a  straight- line  manner  with  a  minimal 
amount  of  branching  for  recurrent  processes.   The  original 
program  processed  data  for  a  least-squares  fitting  of  a 
Gaussian  and  calculated  the  lognormal  variance  (SIGMA)  and 
the  refractive  index  structure  constant  (CJ  but  did  not 
compute  the  average  detector  voltage  or  an  extinction  coef- 
ficient.  The  program  lacked  the  versatility  to  accept  data 
from  available  analyzers  by  manual  input. 

The  Hewlett-Packard  9  810A  calculator  (with  option  003 
installed)  has  a  total  of  2  036  locations  in  memory  available 
for  program  steps  of  which  1847  were  used  in  the  original 
program.   Since  the  additional  desired  processing  could  not 
possibly  be  accomplished  in  the  steps  remaining,  a  major 
program  overhaul  was  required.   The  first  measure  taken  was 
to  delete  or  shorten  as  many  labels  and  operator  instructions 
as  possible;  approximately  2  50  program  steps  were  excised 
in  this  manner  with  no  functional  loss.   More  steps  were 
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saved  by  recording  implemental  constants  on  a  magnetic  card 
and  entering  them  directly  into  the  proper  memory  locations 
immediately  after  the  program  was  loaded. 

The  program  now  had  enough  room  for  additional  processing 
In  previous  analysis  [Ref.  3,  et.al.)  it  was  assumed  that 
the  natural  log  of  the  average  detector  voltage  (beam  inten- 
sity) was  approximately  equal  to  the  average  natural  log  of 
detector  voltage.   The  validity  of  this  assumption,  initially 
made  when  computing  the  log-intensity  variance  with  photo- 
graphs, is  crucial  in  extinction  measurement  (which  depends 
on  accurate  determination  of  the  average  detector  voltage) . 
The  immediate  concern  was  the  possibility  of  a  functional 
relationship  between  average  detector  voltage  (V)  and  the 
exponential  of  the  average  log  detector  voltage  (exp (In (V) ) ) 


dependent  on  the  average  value  (ln(V))  and  the  standard 
deviation  (a)  of  the  lognormal  distribution. 

A  computer  program  was  developed  which  analyzed  the 
variation  of  the  ratio  of  the  two  quantities  with  increments 


of  either  ln(V)  or  a.   The  method  was  to  assume  an  ln(V) 
value  for  a  Gaussian  distribution  and  compute  the  average 
voltage.   The  probability  density  function  of  In (V) was 
assumed,  therefore, to  be  of  the  form: 


_  r     (X   -    X)     , 

Yin  =  c  expt     2   1 


where: 
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x  =  ln(V) 
x  =  In  (v) 

Let  y (V)  be  the  probability  density  function  of  V: 

b  exp(b) 

/      y.    dx  =  /  y(V)    dV 

a  ln  exp(a) 


dx  =  d{ln(V) }   =  ^ 


a  dV       yln    ,„ 

•  •     ym  **  -  Ym  —  =  —  dv 


oo  ln(<»)    yn  » 

/      y,    dx  =  /  -±&  dV  -      /      y(V)dV 

ln  liiO)  V  0 


Since  the  number  of  total  counts  accumulated  over  a  range 
of  detector  voltage  while  sampling  will  be  identical  whether 
the  display  distribution  is  linear  or  logarithmic: 


yln 
y (V) dV  =  ylndX  =  -±£  dV 


Yin  r  ( v_v\ ^ 


...  y(v,  =  *L   =  §  expt-ii^-]  =  §  Bp|.'I'W^1 

2a  2a 


Applying  statistics: 


-  m   /Vy(V)dV 
/y(V)dV 
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For  a  normal  distribution  a  dispersion  index  of  ± 3a  gives 
a  confidence  limit  of  99.73%  [Ref.  7,  et.al.)  and  these 
limits  were  imposed  in  a  program  applying  Weddle  *  s  rule 
[Ref.  8]  to  accomplish  the  integration.   The  ratio  of  V 


to  exp(ln(V))  was  found  to  be  independent  of  ln(V)  but  an 
increasing  function  of  a.   With  points  generated  by  the  pro- 
gram (for  a   values  between  0.1  and  1.2  at  increments  of  0.02) 
a  fourth-degree  polynomial  was  found  (using  least-squares 
fitting)  to  generate  average  detector  voltage  values  from 
average  log-intensity.   The  polynomial  coefficients  were 
entered  on  the  magnetic  card  carrying  program  constants 
and  used  to  find  the  average  detector  voltage  after  deter- 


mining ln(V)  (from  the  coefficients  of  the  least-squares 
polynomial  fitting  data  to  a  lognormal  distribution) .   An 
easier  method  was  subsequently  employed  which  can  be 
justified  by  proceeding  further  with  the  preceding  derivation, 

==  _  /Vy(V)dV 


V  = 


/y  (V)dV 


but 


V  =  exp(ln(V))  =  exp(x) 


and 


y(V)dV  =  yln  dx 
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/exp(x)ylndx   fC   exP(x)  exp[-i^-]dx 


_  2 

/Y-i    dx                       ._,             r     (x-x)     ,     , 
Jjrln  /C   exp[— j1 — ]    dx 

2<j 


After  combining  exponentials  and  completing  the  square  in 

—    2 

the  numberator,  and  defining  x1  =  x  -  a  : 

2  /         M2 

C  exp(x+^-)    /exp[-   x~x        ]    dx 

V  =  -  2a 


C   /exp[-(x~x)     ]    dx 
2oZ 


The  integrands  in  numerator  and  denominator  are  Gaussian 
distributions  with  the  same  standard  deviations  and, 
therefore,  the  same  areas. 


.\  V  =  exp(x  +  %")  =  exp(ln(V))  exp(~) 


This  derived  result  can  be  found  in  meticulous  statistical 


literature  [Ref.  9].   Since  In (V)  and  a  were  already  avail- 
able in  the  program  the  necessary  calculation  was  easily 
included  and  required  significantly  less  memory  than  the 
fourth-degree  polynomial  conversion.   Incidentally,  it  was 
found  that  the  polynomial  method  was  accurate  within  2% 
for  standard  deviation  values  between  0.1  and  0.9! 

With  an  assumed  extinction  constant  the  calculation  of 
an  extinction  coefficient  is  easily  made,  i.e.: 

ln(-%) 

vir 

n  =  5 — 
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where 


u  =  extinction  coefficient 

2 

C  =  extinction  constant  (V  R   ) 

o  o 

V  =  average  detector  voltage 
R  =  range 

This  computation  was  added  to  the  program  with  the  extinction 
constant  becoming  one  of  the  operator  input  parameters. 

The  next  program  modification  was  the  addition  of  a 
section  which  could  process  linear  data  to  find  the  average 
detector  voltage  and  extinction  coefficient.   The  main 
purpose  was  to  correlate  results  with  those  from  the  log- 
arithmic section  but  this  was  not  done  to  any  great  extent. 
The  average  channel  number  is  determined  by  dividing  the 
summation  of  channel  number  times  channel  counts  by  the 
summation  of  counts.   The  average  channel  number  is  then 
modified  in  a  linear  scaling  process  to  give  the  average 
voltage.   The  extinction  coefficient  was  computed  in  exactly 
the  same  way  as  in  the  logarithmic  section. 

After  the  linear  portion  was  added  the  program  almost 
filled  available  memory  and  was  found  to  execute  more  slowly 
than  before  —  especially  plotting  data.   The  program  showed 
minimal  utilization  of  subroutines  and,  consequently,  the 
program  wasted  memory.   The  subroutines  used  were  scattered 
throughout  the  body  of  the  program  which  meant  wasted  time 
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spent  searching  after  branching  statements.   The  computer 
still  took  data  input  only  from  the  PIP-400  through  the 
interface.   A  Nuclear  Data  128  channel  analyzer  which  inter- 
faced with  a  teletype  was  available  as  a  back-up  but  the 
program  had  to  be  modified  to  accept  manual  input  and  a 
different  channel  capacity.   With  all  the  discrepancies 
a  major  reorganization  was  needed. 

The  HP9810A  calculator  branches  to  subroutines  by  going 
to  the  program  beginning  and  searching  through  until  finding 
the  called  label.   This  requires  an  inordinate  amount  of 
time  to  complete  "do-loops"  and  subroutines.   This  problem 
was  solved  by  making  processes  which  were  reiterative  or 
used  by  more  than  one  section  of  the  program  into  subroutines 
and  placing  them  in  front  of  the  main  body.   "Flag"  state- 
ments were  used  to  allow  automatic  or  manual  data  input  from 
the  PIP-400  or  Nuclear  Data  analyzer  respectively. 

The  original  program  used  a  large  amount  of  memory  in  a 
calculation  to  filter  out  bias  introduced  in  the  scintillo- 
meter.  When  the  single-detector  system  for  a  chopped  laser 
source  was  first  put  together  the  demodulator  would  output 
a  few  tenths  of  a  volt  with  no  input  present  thus  biasing 
the  sample  distribution.   This  offset  had  to  be  taken  into 
account  when  processing  data  because  of  its  effect  when 
using  the  logarithmic  scale.   In  subsequent  demodulators 
the  offset  problem  was  solved  and  the  calculation  became 
extraneous.   When  the  process  was  removed  from  the  program 
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the  only  discernable  difference  was  an  increased  speed  in 
plotting  data. 

The  net  result  of  all  the  program  modifications  pre- 
viously discussed  was  a  drop  in  the  average  processing 
time  of  an  experimental  distribution  from  over  ten  minutes 
to  less  than  six  minutes. 

B.       METEOROLOGY    PROGRAMS 

Just  prior  to  the  most  recent  R/V  Acania  experiments 
(May,  1976)  Professor  Schacher  (associated  with  the  atmos- 
pheric turbulence  group)  requested  a  computer  program  to  fit 
a  Gaussian  curve  to  data  manually  input  from  the  Nuclear 
Data  pulse  height  analyzer.   No  explanation  of  the  relevant 
theory  was  given  or  required  to  distill  the  necessary 
processing  from  the  existing  scintillation/extinction  program. 

Upon  analyzing  the  processed  data  it  was  found  that  the 
Gaussian  fit  was  a  poor  one.   Prof.  Schacher  suggested  the 
resultant  data  distributions  appeared  more  Lorentzian  and 
requested  another  program  fitting  data  to  that  distribution. 

The  Lorentzian  distribution  is  of  the  form: 


y  = 


2       —  2       2     —     —2     2 

rz  +  (x-x)     x^  -  2xx  +  xz  +  r 


i  « 4  x2  _  2x   +  x2  +   r2 

y   A       A         A 


This  can  again  be  fitted  by  least-squares  method  to  a 
quadratic  and  the  coefficients  evaluated  to  get  the  average 
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voltage  (after  scaling  x,  the  average  channel  number)  and  r 
(after  scaling) .   The  weighting  scheme  differed  from  that 
used  for  the  Gaussian.   The  uncertainty  is  (again  assuming 
the  number  of  counts  in  each  channel  is  Poisson-distributed) 


,,1,      -2,    -dy   -/y     -3/2 
d(-)  =-y   dy  =  — X  =  _yL  =  -y 

2  y      y 


Again  the  data  was  weighted  inversely  with  its  uncertainty 
(multiplied  by  the  number  of  counts  in  the  channel  to  the 
3/2  power) .   The  least-squares  fit  attempts  to  minimize: 


1       2 

Ego.  ( —  -  ax.   -  bx.  -  c) 
l  yi     i      l 


where 


3/2 

go.  =  y . 


y.  =  number  of  counts  in  channel  i 

2  l 


x.  =  channel  number  i 

l 


1 
a   '  A 


.     -2x 

b  _  _ 

^2+  r2 

C  =  r 


Since  all  the  machinery  of  the  least-squares  fitting  of 
a  quadratic  applying  Cramer's  rule  already  existed,  it  was 
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a  relatively  simple  matter  to  modify  the  existing  Gaussian 
program.   The  comparative  results  definitely  point  to  a 
Lorentzian  distribution,  which  was  not  predicted  by  theory, 
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A  SIMPLIFIED  FLOW  CHART 
OF  THE  NEW  PROGRAM 


YES 


ENTER  CALIBRATION  CHANNELS, 
CALIBRATION  VOLTAGES 


CONSTRUCT  SCALE  FACTOR 


ENTER  EXTINCTION  CONSTANT, 
WAVELENGTH,  RANGE,  PAR  GAIN, 
MAXIMUM  COUNT,  TIME 


LOGARITHMIC 


PLOT  AND  LABEL  AXES 


ENTER  DATA  AND  PLOT.   CON- 
STRUCT SUMMATIONS  FOR  LEAST- 
SQUARES  FIT. 


CONSTRUCT  COEFFICIENTS  OF 
LEAST-SQUARES  QUADRATIC 


CONSTRUCT  LOGARITHMIC 
VOLTAGE  VARIANCE 


PLOTTER  PRINT;  DATE,  TIME, 
LMBD  (WAVELENGTH),  RNGE 

(range),  $  SQ  (variance), 
$  (standard  deviation) 


CONSTRUCT  AND  PLOTTER  PRINT 

CNSQ  (CN  SQUARED),  CN  (c„) 
N N 


CONSTRUCT  AND  PLOTTER  PRINT 
AVE  V  (AVERAGE  VOLTAGE) 


CALIBRATION  POINTS 
PRINTED  OUT 


LINEAR 


ENTER  DATA,  CONSTRUCT  SUMMATIONS 
FOR  AVERAGE  VOLTAGE 


CONSTRUCT  AND  PRINT  AVE  V 
(AVERAGE  voltage) 


CONSTRUCT  AND  PRINT  EXT  MU 
(EXTINCTION  coefficient) 


|  END^- 


CONSTRUCT  AND  PLOTTER  PRINT 
EXT  MU  (EXTINCTION  COEFFICIENT) 
PLOTTER  PRINT  GAIN  (PAR  GAIN) 


PLOTTER  PLOT  COMPUTED  GAUSSIAN 
CURVE  AND  CALIBRATION  MARKERS 


I  end\- 
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00103^+ 

000000 
000000 
000000 
006291 
000000 
000000 
000000 
000000 
006638 
000000 
000000 
000000 
000000 
000000 
000000 


♦000000 
♦000000 
♦000000 
♦000000 
♦008860 
♦000000 
♦000000 
♦000000 
♦000000 
♦003506 
♦000000 
♦000000 
♦000000 
♦000000 
♦000000 
♦000000 


♦000001 
♦000000 
♦000000 
♦000000 

♦005971 

♦000000 

♦000000 

♦000000 

♦000000 

♦000000 

♦000000 

♦000000 

♦000000 

♦000000. 

♦000000 

♦000000 


♦000000 
♦000000 
♦000000 
♦000000 
♦000000 
♦000000 
♦000001 
♦000000 
♦000000 
♦000000 
♦000000 

♦006799 
♦000000 
♦000000 
♦000000 
♦000000 


♦000000 
♦000000 
♦000000 
♦000000 
♦000000 
♦000000 

♦009367, 
♦000000 
♦000000 
♦000000 
♦000000 
♦002922 
♦000000 
♦000000 
♦000000 
♦000000 


♦000001 
♦000000 
♦000000 
♦000000 
♦000000 
♦000000 
♦016270 
♦000000 
♦000000 
♦000000 
♦000000 
♦000000 
♦000000 
♦000121 
♦000000 
♦000000 


♦000000 
♦000000 
♦000000 
♦000000 
♦000000 
♦000000 
♦000002 
♦000000 
♦000000 
♦000000 
♦000000 
♦000000 
♦000000 
♦007691 
♦000000 
♦000000 


♦000000 
♦000000 
♦000000 
♦000000 
♦000000 
♦000001 
♦000000 
♦000000 
♦000007 
♦000000 
♦000000 
♦000000 
♦000000 

♦003369 
♦000000 
♦000000 


An  example  of  the  teletype  output  of  the  voltage  calibration 
channel  numbers  used  in  scaling  (reproeduced  from  1  March  1976) 
Notice  the  extraneous  counts  collected  in  channel  1.   The 
points  would  be  entered  as  follows: 

Voltage 
0.25 
0.50 
1.00 
2.00 
4.00 


Channel  # 

Change  in  # 

32.97 

52.63 

19.60 

72.35 

19.72 

91.30 

18.95 

110.30 

19.00 

PICTURE  OF  EXPERIMENTAL  DISTRIBUTION  AS  DISPLAYED  ON  THE 
PIP-400  PULSE  HEIGHT  ANALYZER 
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SET  FLAG  FOB 
GAL  PRINTOUT 

11701.073 

3117.074 

7040.103 

10768.104 

10883.134 

10939.164 

8756.194 

2194.195 
ENTER  1ST  CaL 
ENTER  2ND.  GAL 
ENTER  3RD  GAL 
ENTER  VI  IN  I 
AND  U3  IN  X 
ENTER.  VOLTAGE 
OFFSET- 
ENTER,  WAVELENGTH 
IN  MICRONS 
ENTER.  RANGE 
IN  METERS, 
ENTER.  DATE 
ENTER  TIME 
SET.  ANALYZER.  IQ 
PRINT.  OUT.  AND. 
SET.  FLAG  FOR 
DATA  PRINT  OUT 
ENTER  MAX  COUNT- 


DATE 
TIME 

WAVELENGTH 
RANGE 


20776. 

1245. 

0.6328 

900. 


183.092 
291.093 
503.094 

569.095 
645.096 

640.097 
579.098 
758.099 
714.100 
767.101 
849.102 
(etc) 


MONTH? 

4. 

DAY? 

30. 

YEAR? 

76. 

FLAG  -  CAL  PTS 
1ST  CAL? 

61.0 
2ND  CAL? 

117.0 
3RD  CAL? 

173.0 
VI -Y  AND  V3-X 

0.50 
8.00 
EXTINCT  CONST? 

15300.00 
WAVELENGTH 
(MICRONS)? 

0.6328 
RANGE  (METERS)? 

1625. 
GAIN? 

200. 
MAX  COUNT? 

8000. 
TIME? 

1311. 
FLAG  -  LINEAR 
FLAG  -  LOAD  DATA 
PIP  ON  PRINT? 
END, 


MONTH? 

DAY? 

YEAR? 


4. 
28. 
76. 


1ST  CAL? 
2ND.  CAL? 
3RD  CAL? 
VI -Y  AND  V3-X 
MAX  COUNT? 

80000.00 
TIME? 

1555.00 
LOAD  DATA 
END 


Reproductions  of  computer  printer 
processing:  original  program  with 
(left),  present  program  (upper  rt. 


tapes  made  during  data 
cal  and  data  printout 
) ,  meteorology  program 
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APPENDIX  A 


The  following  lists  the  connections  made  in  mating  the 
output  plug  of  the  PIP-400  to  the  input  plug  of  the  interface 
for  the  HP9810A: 

CALCULATOR 
LEAD  FUNCTION         PIP  LEAD  #        LEAD  #  COLOR 


DATA 


10      - 


1CT    - 


ADDRESS 


io- 


10' 


DATA 


10' 


bit  1 
"  2 
"  4 
"       8 

1  bit 

2  " 
4  " 
8      " 


10      -   1 
2 


-  1 
2 
4 
8 

-  1 
2 
4 
8 


-    1 
2 

4 

2         8 
10      -    1 

2 

4 

1         8 

lcr  -  i 

2 

4 

io0-? 

2 

4 
8 
Print  Command 
Hold  Command 
Ground 

+  6  volts 


bit 


bit 


bit 
ti 

ti 

it 

bit 


bit 

it 
it 

bit 
ii 

it 

ti 

bit 
n 

ti 

■I 


1  6 

2  7 

3  8 

4  9 

5  10 

6  11 

7  12 

8  13 
9-12  N/C 

13  30 

14  31 
15,16  N/C 

17  34 

18  35 

19  36 

20  37 

21  38 

22  39 

23  40 

24  41 

25  14 

26  15 

27  16 

28  17 

29  18 

30  19 

31  20 

32  21 

33  22 

34  23 

35  24 

36  25 

37  26 

38  27 

39  28 

40  29 

41  inverted   4  9 

42  ■  47 
43 

44-4 
49 
50 


red 
white/black 
white/green 

blue 
white/blue 

green 
white/gray 

yellow 


white 
white 


white 


white 


white/violet 
white/black (2) 
white/yellow/black 
white/blue/black 
white/red/black 
white/orange/black 

black 
white/gray/black 

white/yellow 
white/brown/black 

white/orange 
white/green/black 
white/violet/black 
white/red 
white/brown 
white/black 
white 


chassis  and  power  supply  white 
8  N/C 

hex  inverter      white 
N/C 


N/C  designates  leads  which  are  not  connected 
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APPENDIX  B 

The  following  is  a  listing  of  the  data  storage  registers 
utilized  by  the  program: 

USE 
computed  average  logarithm  of  detector  voltage 


computed  summations  for  least-squares  fit;  001, 
002  and  003  are  used  subsequently  to  store  the 
coefficients  of  the  least-squares  polynomial 


9999 

determinants  (formed  from  summations)  which  give 
coefficients  of  the  least-squares  polynomial 

computed  standard  deviation  of  least-squares 
Gaussian  fitted  to  data 
first  calibration  channel  number 
second 
third 

computed  variance  of  least-squares  Gaussian, 
later  used  again  in  C.T  calculation 
time 

wavelength  in  meters 
range  in  meters 
wavelength  in  microns 
used  in  CN  calculation 
month 
day 
year 

5/max  count* 
gain 

extinction  constant 
computed  average  detector  voltage 
computed  extinction  constant 
data  point 

data  point  weighting 

channel  number-to-voltage  scaling  constant 
do-loop  accumulators  for  the  linear  calculation 
of  average  detector  voltage 
first  calibration  voltage 
third 

varying  pen  coordinate  value  for  axis  labelling 
incremental  pen  coordinace  values  for  axis  labelling 
label  values  in  axis  labelling 
236  (range  of  x-axis) 
11.8  (range  of  y-axis) 
alpha  scale  factor  (0.1,  FMT  11) 
computer-scaled  x-coordinate 
"         y-coordinate 
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REGISTER 

000 

001 

002 

003 

004 

005 

006 

007 

008 

009** 

010 

011 

012 

013 

014 

021* 

022* 

023* 

025 

028* 

030 

031* 

041* 

050 

051* 

052* 

053* 

054 

055* 

056* 

057 

058 

060 

061 

065 

068 

069 

078* 

079* 

095 

097 

098 

099** 

100** 

102** 

103 

104 

REGISTER 

USE 

105** 

-0.9 

(Y  .  ) 

106** 

10.9 

jYminj 

107** 

-18 

<X^n> 

108** 

218 

max 

(followed  by  FMT  13) 
(followed  by  FMT  12) 
a  &  b  used  many  times  through  out  the  program 

*   marks  those  registers  whose  values  are  input  by  the  operator 
**  marks  those  registers  with  values  stored  on  a  magnetic  card 


The  following  is  a  list  of  label  titles  in  the  order 
they  appear  in  the  program: 

LABEL  FUNCTION 

takes  manual  or  interface  input;  puts  data  in 

060,  channel  in  a,  5/max  count  in  b 

puts  channel  number  in  y  and  upper  data  channel 

limit  in  x 

puts  channel  number  in  y  and  upper  channel 

limit  in  x 

plots  data  points  and  constructs  summations 

necessary  for  least-squares  fitting 

plots  the  computed  Gaussian  and  calibration  points 

plots  marker  at  calibration  point 

constructs  a  special  floating  point  form  of 

a  number  to  be  plotted 

labels  plotted  x  &  y  axes 

computes  the  extinction  coefficient 

prints  operator  instructions: 

FLAG  -  LOAD  DATA;  PIP  ON  PRINT? 
lists  calibration  points 

preliminary  program  set-up  (month,  day,  year, 
calibration  points,  extinction  constant) 
completes  program  set-up  (wavelength,  range, 
gain,  max  count,  time);  branches  to  label  N 
for  linear;  plots  and  labels  x  &  y  axes;  plots 
data;  calculates  desired  constants;  plots  con- 
stants, computed  Gaussian,  and  calibration 
points;  returns  to  label  M 
N(J,A,B,I,      linear  portion  of  program;  computes  average 
M,8,9)       detector  voltage  and  extinction  constant, 
and  prints  them,  returns  to  label  M 

*   marks  labels  which  are  subroutines 

()  hold  subroutines  and  labels  branched  to 


A*  (UP 

,DN) 

B*(l, 

2) 

C*(3, 

4) 

D*(A, 

B,C,5) 

E*(F) 

F* 

G* 

H*(6, 

I* 

J* 

7) 

K* 

L(K) 

M(N,H 
G,I 

,J,D, 
,E) 
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APPENDIX  C 

The  following  is  a  complete  listing  of  the  program  as 

printed  out  by  the  calculator; 

0000  GTO  44  0050  LBL  51  0100    +  33 

0001  LBL  51  0051  DN  25  0101    DN  25 

0002  L  72  0052  LBL  51  0102  CNT  47 

0003  LBL  51  0053  UP  27  0103  S/R  77 

004  A  62  0054  FMT  42  0104  LBL  51 

0005  IFG  43  0055  1  01  0105    B  66 

0006  SFL  54  0056  8  10  0106    a  13 

0007  GTO  44  0057  0  00  0107    UP  27 

0008  LBL  .   51  0058  1  01  0108  IFG  43 

0009  UP  27  0059  XTO  23  0109  SFL  54 

0010  FMT  42  0060  +  33  0110  GTO  44 

0011  1  01  0061  0  13  0111  LBL  51 

0012  8  10  0062  0  13  0112    1  01 

0013  3  03  0063  UP  27  0113  CNT  47 

0014  FMT  42  0064  CLX  37  0114    1  01 

0015  3  03  0065  STP  41  0115    9  11 

0016  3  03  0066  XTO  23  0116    5  05 

0017  .  21  0067  6  06  0117  GTO  44 

0018  UP  27  0068  0  00  0118  LBL  51 

0019  CLX  37  0069  FMT  42  0119    2  02 

0020  X=Y  50  0070  1  01  0120  LBL  51 

0021  GTO  44  0071  8  10  0121    1  01 

0022  LBL  51  0072  3  03  0122  CNT  47 

0023  A  62  0073  LBL  51  0123    1  01 

0024  CNT  47  0074  DN  25  0124    2  02 

0025  EEX  26  0075  XFR  67  0125    3  03 

0026  CHS  32  0076  5  05  0126  LBL  51 

0027  3  03  0077  4  04  0127    2  02 

0028  X  36  0078  XTO  23  0128  S/R  77 

0029  DN  25  0079  b  14  0129  LBL  51 

0030  UP  27  0080  1/X  17  0130    C  61 

0031  INT  64  0081  UP  27  0131    a  13 

0032  -  34  0082  1  01  0132    UP  27 

0033  XTO  23  0083  0  00  0133  IFG  43 

0034  6  06  0084  DIV  35  0134  SFL  54 

0035  0  00  0085  XFR  67  0135  FTO  44 

0036  EEX  26  0086    6  06  0136  LBL  51 

0037  3  03  0087    0  00  0137    3  03 

0038  X  36  0088  X<Y  52  0138    1  01 

0039  2  02  0089  CNT  47  0139    9  11 

0040  0  00  0090  S/R  77  0140    9  11 

0041  0  00  0091  CNT  47  0141  GTO  44 

0042  X<Y  52  0092  CNT  47  0142  LBL  51 

0043  CNT  47  0093  UP  27  0143    4  04 

0044  -  34  0094    a  13  0144  LBL  51 

0045  CNT  47  0095  UP  27  0145    3  03 

0046  CNT  47  0096  EEX  26  0146    1  01 

0047  YTO  40  0097    3  03  0147    2  02 

0048  0  13  0098  DIV  35  0148    8  10 

0049  GTO  44  0099    DN  25  0149  LBL  51 
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0150 

4 

04 

0200 

GTO 

44 

0250 

4 

04 

0151 

S/R 

77 

0201 

LBL 

51 

0251 

XFR 

67 

0152 

LBL 

51 

0202 

5 

05 

0252 

X 

36 

0153 

D 

63 

0203 

CNT 

47 

0253 

a 

13 

0154 

FTO 

44 

0204 

XFR 

67 

0254 

XTO 

23 

0155 

S/R 

77 

0205 

6 

06 

0255 

+ 

33 

0156 

LBL 

51 

0206 

0 

00 

0256 

5 

05 

0157 

A 

62 

0207 

I 

65 

0257 

LBL 

51 

0158 

XFR 

67 

0208 

XFR 

67 

0258 

5 

05 

0159 

6 

06 

0209 

X 

36 

0259 

GTO 

44 

0160 

0 

00 

0210 

6 

06 

0260 

S/R 

77 

0161 

UP 

27 

0211 

1 

01 

0261 

LBL 

51 

0162 

S~ 

76 

0212 

XTO 

23 

0262 

C 

61 

0163 

XTO 

23 

0213 

+ 

33 

0263 

X>Y 

53 

0164 

6 

06 

0214 

6 

06 

0264 

GTO 

44 

0165 

1 

01 

0215 

XFR 

67 

0265 

LBL 

51 

0166 

b 

14 

0216 

X 

36 

0266 

D 

63 

0167 

X 

36 

0217 

a 

13 

0267 

CNT 

47 

0168 

a 

13 

0218 

XTO 

23 

0268 

S/R 

77 

0169 

FMT 

42 

0219 

+ 

33 

0269 

LBL 

51 

0170 

1 

01 

0220 

7 

07 

0270 

E 

60 

0171 

UP 

27 

0221 

XFR 

67 

0271 

a 

13 

0172 

FMT 

42 

0222 

X 

36 

0272 

UP 

27 

0173 

1 

01 

0223 

a 

13 

0273 

XSQ 

12 

0174 

4 

04 

0224 

XTO 

23 

0274 

XFR 

67 

0175 

CNT 

47 

0225 

+ 

33 

0275 

X 

36 

0176 

5 

05 

0226 

8 

10 

0276 

1 

01 

0177 

UP 

27 

0227 

XFR 

67 

0277 

YE 

24 

0178 

a 

13 

0228 

6 

06 

0278 

X 

36 

0179 

X<Y 

52 

0229 

1 

01 

0279 

2 

02 

0180 

GTO 

44 

0230 

XTO 

23 

0280 

XFR 

67 

0181 

LBL 

51 

0231 

+ 

33 

0281 

+ 

33 

0182 

5 

05 

0232 

1 

01 

0282 

3 

03 

0183 

CNT 

47 

0233 

XFR 

67 

0283 

+ 

33 

0184 

GTO 

44 

0234 

X 

36 

0284 

DN 

25 

0185 

S/R 

77 

0235 

a 

13 

0285 

J 

75 

0186 

LBL 

51 

0236 

XTO 

23 

0286 

XFR 

67 

0187 

B 

66 

0237 

+ 

33 

0287 

X 

36 

0188 

X<Y 

52 

0238 

2 

02 

0288 

b 

14 

0189 

GTO 

44 

0239 

XFR 

67 

0289 

UP 

27 

0190 

LBL 

51 

0240 

X 

36 

0290 

a 

13 

0191 

5 

05 

0241 

a 

13 

0291 

FMT 

42 

0192 

CNT 

47 

0242 

XTO 

23 

0292 

1 

01 

0193 

XFR 

67 

0243 

+ 

33 

0293 

DN 

25 

0194 

6 

06 

0244 

3 

03 

0294 

1 

01 

0195 

0 

00 

0245 

XFR 

67 

0295 

XTO 

23 

0196 

UP 

27 

0246 

X 

36 

0296 

+ 

33 

0197 

b 

14 

0247 

a 

13 

0297 

a 

13 

0198 

l/x 

17 

024  8 

XTO 

23 

0298 

a 

13 

0199 

X>Y 

53 

0249 

+ 

3 

0299 

UP 

27 

59 


0300 

2 

02 

0350 

- 

34 

0400 

9 

11 

0301 

0 

00 

0351 

UP 

27 

0401 

8 

10 

0302 

0 

00 

0352 

1 

01 

0402 

XTO 

23 

0303 

X<Y 

53 

0353 

- 

34 

0403 

1 

01 

0304 

GTO 

44 

0354 

RUP 

22 

04C4 

0 

00 

0305 

LBL 

51 

>  0355 

+ 

33 

0405 

4 

04 

0306 

E 

60 

0356 

DN 

25 

0406 

XFR 

67 

0307 

CNT 

47 

0357 

K 

55 

0407 

9 

11 

0308 

XFR 

67 

0358 

5 

05 

0408 

5 

05 

0309 

2 

02 

0359 

FMT 

42 

0409 

FMT 

42 

0310 

1 

01 

0360 

1 

01 

0410 

1 

01 

0311 

UP 

27 

0361 

8 

10 

0411 

PNT 

45 

0312 

F 

16 

0362 

3 

03 

0412 

XFR 

67 

0313 

XFR 

67 

0363 

YTO 

40 

0314 

9 

11 

0314 

2 

02 

0364 

5 

05 

0414 

7 

07 

0315 

2 

02 

0365 

0 

00 

0415 

XTO 

23 

0316 

UP 

27 

0366 

FMT 

42 

0416 

+ 

33 

0317 

F 

16 

0367 

1 

01 

0417 

1 

01 

0318 

XFR 

67 

0368 

PNT 

45 

0418 

0 

00 

0319 

2 

02 

0369 

FMT 

42 

0419 

3 

03 

0320 

3 

03 

0370 

1 

01 

0420 

1 

01 

0321 

UP 

27 

0371 

FMT 

42 

0421 

XTO 

23 

0322 

F 

16 

0372 

CNT 

47 

0422 

+ 

33 

0323 

S/R 

77 

0373 

FMT 

42 

0423 

9 

11 

0324 

LBL 

51 

0374 

FMT 

42 

0424 

5 

05 

0325 

F 

16 

0375 

1 

01 

0425 

GTO 

44 

0326 

CLX 

37 

0376 

8 

10 

0426 

LBL 

51 

0327 

XEY 

30 

0377 

0 

00 

0427 

H 

74 

0328 

FMT 

42 

0378 

XFR 

67 

0428 

LBL 

51 

0329 

1 

01 

0379 

5 

05 

0429 

6 

06 

0330 

UP 

27 

0380 

0 

00 

0430 

CLR 

20 

0331 

FMT 

42 

0381 

FMT 

42 

0431 

FMT 

42 

0332 

DN 

25 

0382 

1 

01 

0432 

1 

01 

0333 

UP 

27 

0383 

PNT 

45 

0433 

UP 

27 

0334 

1 

01 

0384 

S/R 

77 

0434 

1 

01 

0335 

5 

05 

0385 

LBL 

51 

0435 

9 

11 

0336 

0 

00 

0386 

H 

74 

0436 

0 

00 

0337 

0 

00 

0387 

XFR 

67 

0437 

XTO 

23 

0338 

XEY 

30 

0388 

1 

01 

0438 

- 

34 

0339 

FMT 

42 

0389 

0 

00 

0439 

1 

01 

0340 

DN 

25 

0390 

3 

03 

0440 

0 

00 

0341 

FMT 

42 

0391 

UP 

27 

0441 

3 

03 

0342 

UP 

27 

0392 

XFR 

67 

0442 

XFR 

67 

0343 

S/R 

77 

0393 

9 

11 

0443 

1 

01 

0344 

LBL 

51 

0394 

X<Y 

52 

0444 

0 

00 

0345 

G 

15 

0395 

GTO 

44 

0445 

4 

04 

0346 

K 

55 

0396 

LBL 

51 

0446 

XTO 

23 

0347 

4 

04 

0397 

6 

06 

0447 

0 

00 

0348 

UP 

27 

0398 

CNT 

47 

0448 

9 

11 

0349 

INT 

64 

0399 

XFR 

67 

0449 

8 

10 

60 


0450 

2 

02 

0500 

XTO 

23 

0550 

0 

71 

0451 

0 

00 

0501 

+ 

33 

0551 

A 

62 

0452 

UP 

27 

0502 

9 

11 

0552 

D 

63 

0453 

XFR 

67 

0503 

5 

05 

0553 

CNT 

47 

0454 

9 

11 

0504 

GTO 

44 

0554 

D 

63 

0455 

9 

11 

0505 

LBL 

51 

0555 

A 

62 

0456 

DIV 

35 

0605 

7 

07 

0556 

XTO 

23 

0457 

XFR 

67 

0507 

S/R 

77 

0557 

A 

62 

0458 

9 

11 

0508 

LBL 

51 

0558 

FMT 

42 

0459 

X 

36 

0509 

I 

65 

0559 

STP 

41 

0460 

YTO 

40 

0510 

XFR 

67 

0560 

IFG 

43 

0461 

9 

11 

0511 

5 

05 

0561 

CNT 

47 

0462 

7 

07 

0512 

7 

07 

0562 

SFL 

54 

0463 

CLX 

37 

0513 

XFR 

67 

0563 

S/R 

77 

0464 

XTO 

23 

0514 

DIV 

35 

0564 

CNT 

47 

0465 

9 

11 

0515 

5 

05 

0565 

FMT 

42 

0466 

5 

05 

0516 

5 

05 

0566 

FMT 

42 

0467 

LBL 

51 

0517 

UP 

27 

0567 

IT 

56 

0468 

7 

07 

0518 

XFR 

67 

0568 

I 

65 

0469 

XFR 

67 

0519 

3 

03 

0569 

IT 

56 

0470 

9 

11 

0520 

1 

01 

0570 

CNT 

47 

0471 

8 

10 

0521 

XSQ 

12 

0571 

O 

71 

0472 

XTO 

23 

0522 

X 

36 

0572 

N 

73 

0473 

1 

01 

0523 

DN 

25 

0573 

CNT 

47 

0474 

0 

00 

0524 

l/x 

17 

0574 

TT 

56 

0475 

4 

04 

0525 

XFR 

67 

0575 

a 

13 

0476 

UP 

27 

0526 

X 

36 

0576 

I 

65 

0477 

XFR 

67 

0527 

5 

05 

0577 

N 

73 

0478 

9 

11 

0528 

6 

06 

0578 

XTO 

23 

0479 

X<Y 

52 

0529 

1 

65 

0579 

IFG 

43 

0480 

CNT 

47 

0530 

XFR 

67 

0580 

FMT 

42 

0481 

S/R 

77 

0531 

DIV 

35 

0581 

STP 

41 

0482 

CNT 

47 

0532 

3 

03 

0582 

S/R 

77 

0483 

CNT 

47 

0533 

1 

01 

0583 

LBL 

51 

0484 

XFR 

67 

0534 

XTO 

23 

0584 

K 

55 

0485 

9 

11 

0535 

5 

05 

0585 

FMT 

42 

0486 

5 

05 

0536 

8 

10 

0586 

1 

01 

0487 

FMT 

42 

0537 

S/R 

77 

0587 

8 

10 

0488 

1 

01 

0538 

LBL 

51 

0588 

3 

03 

0489 

PNT 

45 

0539 

J 

75 

0589 

FMT 

42 

0490 

XFR 

67 

0540 

FMT 

42 

0590 

3 

03 

0491 

9 

11 

0541 

FMT 

42 

0591 

3 

03 

0492 

7 

07 

0542 

F 

16 

0592 

• 

21 

0493 

XTO 

23 

0543 

L 

72 

0593 

UP 

27 

0494 

+ 

33 

0544 

A 

62 

0594 

EEX 

26 

0495 

0 

00 

0545 

G 

15 

0595 

3 

03 

0496 

9 

11 

J  0546 

CNT 

47 

0596 

DIV 

35 

0497 

8 

10 

0547 

EEX 

26 

0597 

UP 

27 

0498 

2 

02 

0548 

CNT 

47 

0598 

DN 

25 

0499 

0 

00 

0549 

L 

72 

0599 

X<Y 

52 

61 


0600 

DN 

25 

0650 

D 

63 

0700 

S/R 

77 

0601 

PNT 

45 

0651 

A 

62 

0701 

LBL 

51 

0602 

CNT 

47 

0652 

XFR 

67 

0702 

K 

55 

0603 

CNT 

47 

0653 

IFG 

43 

0703 

FMT 

42 

0604 

DN 

25 

0654 

FMT 

42 

0704 

1 

01 

0605 

INT 

64 

0655 

XFR 

67 

0705 

8 

10 

0606 

- 

34 

0656 

5 

05 

0706 

1 

01 

0607 

• 

21 

0657 

2 

02 

0707 

FMT 

42 

0608 

2 

02 

0658 

STP 

41 

0708 

FMT 

42 

0609 

X<Y 

52 

0659 

XTO 

23 

0709 

1 

01 

0610 

- 

34 

0660 

5 

05 

0710 

YTO 

40 

0611 

CNT 

47 

0661 

2 

02 

0711 

XTO 

23 

0612 

CNT 

47 

0662 

PNT 

45 

0712 

CNT 

47 

0613 

CNT 

47 

0663 

FMT 

42 

0713 

C 

61 

0614 

• 

21 

0664 

FMT 

42 

0714 

A 

62 

0615 

1 

01 

0665 

XFR 

67 

0715 

L 

72 

0616 

9 

11 

0666 

E 

60 

0716 

IFG 

43 

0617 

9 

11 

0667 

A 

62 

0717 

FMT 

42 

0618 

X<Y 

53 

0668 

a 

13 

0718 

XFR 

67 

0619 

GTO 

44 

0669 

IFG 

43 

0719 

2 

02 

0620 

LBL 

51 

0670 

FMT 

42 

0720 

1 

01 

0621 

K 

55 

0671 

XFR 

67 

0721 

STP 

41 

0622 

CNT 

47 

0672 

5 

05 

0722 

XTO 

23 

0623 

S/R 

77 

0673 

3 

03 

0723 

2 

02 

0624 

LBL 

51 

0674 

STP 

41 

0724 

1 

01 

0625 

L 

72 

0675 

XTO 

23 

0725 

PNT 

45 

0626 

FMT 

42 

0676 

5 

05 

0726 

FMT 

42 

0627 

1 

01 

0677 

3 

03 

0727 

FMT 

42 

0628 

8 

10 

0678 

PNT 

45 

0728 

2 

Q2 

0629 

0 

00 

0679 

CLX 

37 

0729 

N 

73 

0630 

CLR 

20 

0680 

FMT 

42 

0730 

D 

63 

0631 

FMT 

42 

0681 

FMT 

42 

0731 

CNT 

47 

0632 

FMT 

42 

0682 

F 

16 

0732 

C 

61 

0633 

M 

70 

0683 

L 

72 

0733 

A 

62 

0634 

0 

71 

0684 

A 

62 

0734 

L 

72 

0635 

N 

73 

0685 

G 

15 

0735 

IFG 

43 

0636 

XTO 

23 

0686 

CNT 

47 

0736 

FMT 

42 

0637 

H 

74 

0687 

EEX 

26 

0737 

XFR 

67 

0638 

IFG 

43 

0688 

CNT 

47 

0738 

2 

02 

0639 

FMT 

42 

0689 

C 

61 

0739 

2 

02 

0640 

XFR 

67 

0690 

A 

62 

0740 

STP 

41 

0641 

5 

05 

0691 

L 

72 

0741 

XTO 

23 

0642 

1 

01 

0692 

CNT 

47 

0742 

2 

02 

0643 

STP 

41 

0693 

7T 

56 

0743 

2 

02 

0644 

XTO 

23 

0694 

XTO 

23 

0744 

PNT 

45 

0645 

5 

05 

0695 

YTO 

40 

0745 

FMT 

42 

0646 

1 

01 

0696 

FMT 

42 

0746 

FMT 

42 

0647 

PNT 

45 

0697 

STP 

41 

0747 

3 

03 

0648 

FMT 

42 

0698 

IFG 

43 

0748 

a 

13 

0649 

FMT 

42 

0699 

GTO 

44 

0749 

D 

63 

62 


0750  CNT  47  0800  7  07  0850  PNT  45 

0751  C  61  0801  9  11  0851  LBL  51 

0752  A  62  0802  PNT  45  0852  M  70 

0753  L  72  0803  XFR  67  0853  CLR  20 

0754  IFG  43  0804  7  07  0854  FMT  42 

0755  FMT  42  0805  8  10  0855  FMT  42 

0756  XFR  67  0806  UP  27  0856  IND  31 

0757  2  02  0807  XFR  67  0857  A  62 

0758  3  03  0808  7  07  0858  INT  64 

0759  STP  41  0809  9  11  0859  E  60 

0760  XTO  23  0810  I  65  0860  L  72 

0761  2  02  0811  XEY  30  0861  E  60 

0762  3  03  0812  I  65  0862  N  73 

0763  PNT  45  0813  -  34  0863  G  15 

0764  FMT  42  0814  XFR  67  0864  XTO  23 

0765  FMT  42  0815  2  02  0865  H  74 

0766  INT  64  0816  3  03  0866  CLR  20 

0767  1  01  0817  XFR  67  0867  X<Y  52 

0768  EEX  26  0818  -  34  0868  M  70. 

0769  XFR  67  0819  2  02  0869  I  65 

0770  CNT  47  0820  1  01  0870  C  61 

0771  A  62  0821  DIV  35  0871  a  13 

0772  N  73  0822  YTO  40  0872  0  71 

0773  D  63  0823  6  06  0873  N  73 

0774  CNT  47  0824  5  05  0874  YTO  40 

0775  INT  64  0825  CLR  20  0875  PSE  57 

0776  3  03  0826  FMT  42  0876  IFG  43 

0777  EEX  26  0827  FMT  42  0877  FMT  42 

0778  YE  24  0828  E  60  0878  FMT  42 

0779  FMT  42  0829  YE  24  0879  1  01 

0780  FMT  42  0830  XTO  23  0880  8  10 

0781  1  01  0831  1  65  0881  4  04 

0782  8  10  0832  N  73  0882  XFR  67 

0783  2  02  0833  C  61  0883  4  04 

0784  YE  24  0834  XTO  23  0884  1  .01 

0785  7  07  0835  CNT  47  0885  STP  41 

0786  8  10  0836  C  61  0886  XTO  23 

0787  XFR  67  0837  0  71  0887  4  04 

0788  7  07  0838  N  73  0888  1  01 

0789  9  11  0839  YTO  40  0869  PNT  45 

0790  STP  41  0840  XTO  23  0890  UP  27 

0791  YTO  40  0841  IFG  43  0891  EEX  26 

0792  7  07  0842  FMT  42  0892  6  06 

0793  8  10  0843  XFR  67  0893  DIV  35 

0794  XTO  23  0844  5  05  0894  YTO  40 

0795  7  07  0845  6  06  0895  3  03 

0796  9  11  0846  STP  41  0896  0  00 

0797  XEY  30  0847  XTO  23  0897  CLR  20 

0798  PNT  45  0848  5  05  0898  FMT  42 

0799  XFR  67  0849  6  06  0899  FMT  42 
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0900 

a 

13 

0950 

C 

61 

1000 

N 

73 

0901 

A 

62 

0951 

0 

71 

1001 

E 

60 

0902 

N 

73 

0952 

l/x 

17 

1002 

A 

62 

0903 

G 

15 

0953 

N 

73 

1003 

a 

13 

0904 

E 

60 

0954 

XTO 

23 

1004 

FMT 

42 

0905 

CNT 

47 

0955 

IFG 

43 

1005 

STP 

41 

0906 

X<Y 

52 

0956 

FMT 

42 

1006 

IFG 

43 

0907 

M 

70 

0957 

CLR 

20 

1007 

GTO 

44 

0908 

E 

60 

0958 

5 

05 

1008 

LBL 

51 

0909 

XTO 

23 

0959 

EEX 

26 

1009 

N 

73 

0910 

E 

60 

0960 

3 

03 

1010 

CNT 

47 

0911 

a 

13 

0961 

STP 

41 

1011 

CLR 

20 

0912 

YTO 

40 

0962 

UP 

27 

1012 

FMT 

42 

0918 

PSE 

57 

0963 

PNT 

45 

1013 

1 

01 

0914 

IFG 

43 

0964 

DN 

25 

1014 

UP 

27 

0915 

FMT 

42 

0965 

l/x 

17 

1015 

2 

02 

0916 

FMT 

42 

0966 

UP 

27 

1016 

0 

00 

0917 

1 

01 

0967 

5 

05 

1017 

FMT 

42 

0918 

8 

10 

0968 

X 

36 

1018 

1 

01 

0919 

0 

00 

0969 

YTO 

40 

1019 

5 

05 

0920 

XFR 

67 

0970 

5 

05 

1020 

1 

01 

0921 

3 

03 

0971 

4 

04 

1021 

FMT 

42 

0922 

1 

01 

0972 

CLR 

20 

1022 

1 

01 

0923 

STP 

41 

0973 

FMT 

42 

1023 

6 

06 

0924 

XTO 

23 

0974 
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42 

1024 

4 

04 

0925 

3 

03 

0975 

XTO 

23 

1025 

0 

00 

0926 

1 

01 

0976 

I 
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1026 

0 

00 

0927 

PNT 

45 

0977 

M 

70 

1027 

XTO 

23 

0928 

FMT 

42 

0978 

E 

60 

1028 

0 

00 

0929 

FMT 

42 

0979 

IFG 

43 

1029 

9 

11 

0930 

G 

15 

0980 

FMT 

42 

1030 

8 

10 

0931 

A 

62 

0981 

XFR 

67 

1031 

XFR 

67 

0932 

I 

65 

0982 

2 

02 

1032 

9 

11 

0933 

N 

73 

0983 

8 

10 

1033 

XFR 

67 

0934 

IFG 

43 

0984 

STP 

41 

1034 

DIV 

35 

0935 

FMT 

42 

0985 

XTO 

23 

1035 

1 

01 

0936 

XFR 

67 

0986 

2 

02 

1036 

0 

00 

0937 

5 

05 

0987 

8 

10 

1037 

0 

00 

0938 

5 

05 

0988 

PNT 

45 

1038 

XTO 

23 

0939 

STP 

41 

0989 

FMT 

42 

1039 

9 

11 

0940 

XTO 

23 

0990 

FMT 

42 

1040 

7 

07 

0941 

5 

05 

0991 

F 

16 

1041 

FMT 

42 

0942 

5 

05 

0992 

L 

72 

1042 

1 

01 

0943 

PNT 

45 

0993 

A 

62 

1043 

8 

10 

0944 

FMT 

42 

0994 

G 

15 

1044 

0 

00 

0945 

FMT 

42 

0995 

CNT 

47 

1045 

CLX 

37 

0946 

M 

70 

0996 

EEX 

26 

1046 

XTO 

23 

0947 

A 

62 

0997 

CNT 

47 

1047 

9 

11 

0948 

YE 

24 

0998 

L 

72 

1048 

5 

05 

0949 

CNT 

47 

0999 

I 

65 

1049 

GTO 

44 
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1050  S/R  77  1100  +  33  1150  X  36 

1051  LBL  51  1101  +  33  1151  7  07 

1052  H  74  1102  XFR  67  1152  +  33 

1053  CLR  20  1103  3  03  1153  XFR  67 

1054  XTO  23  1104  XSQ  12  1154  3  03 

1055  1  01  1105  XFR  67  1155  XSQ  12 

1056  XTO  23  1106  X  36  1156  XFR  67 

1057  2  02  1107  3  03  1157  X  36 

1058  XTO  23  1108  -  34  1158  6  06 

1059  3  03  1109  XFR  67  1159  -  34 

1060  XTO  23  1110  2  02  1160  XFR  67 

1061  4  04  1111  XSQ  12  1161  2  02 

1062  XTO  23  1112  XFR  67  1162  XSQ  12 

1063  5  05  1113  X  36  1163  XFR  67 

1064  XTO  23  1114  5  05  1164  X  36 

1065  6  06  1115  -  34  1165  8  10 

1066  XTO  23  1116  XFR  67  1166  -  34 

1067  7  07  1117  4  04  1167  XFR  67 

1068  XTO  23  1118  XSQ  12  1168  1  01 

1069  8  10  1119  XFR  67  1169  XFR  67 

1070  CNT  47  1120  X  36  1170  X  36 

1071  GTO  44  1121  1  01  1171  4  04 

1072  S/R  77  1122  -  34  1172  XFR  67 

1073  LBL  51  1123  YTO  40  1173  X  36 

1074  J  75  1124  1  01  1174  7  07 

1075  GTO  44  1125  0  00  1175  -  34 

1076  S/R  77  1126  XFR  67  1176  YTO  40 

1077  LBL  51  1127  8  10  1177  1  01 

1078  D  63  1128  XFR  67  1178  1  01 

1079  FMT  42  1129  X  36  1179  XFR  67 

1080  1  01  1130  3  03  1180  1  01 

1081  8  10  1131  XFR  67  1181  XFR  67 

1082  0  00  1132  X  36  1182  X  36 

1083  XFR  67  1133  1  01  1183  5  05 

1084  5  05  1134  UP  27  1184  XFR  67 

1085  XFR  67  1135  XFR  67  1185  X  36 

1086  X  36  1136  4  04  1186  7  07 

1087  3  03  1137  XFR  67  1187  UP  27 

1088  XFR  67  1138  X  36  1188  XFR  67 

1089  X  36  1139  2  02  1189  8  10 

1090  1  01  1140  XFR  67  1190  XFR  67 

1091  UP  27  1141  X  36  1191  X  36 

1092  XFR  67  1142  6  06  1192  2  02 

1093  4  04  1143  +  33  1193  XFR  67 

1094  XFR  67  1144  XFR  67  1194  X  36 

1095  X  36  1145  2  02  1195  3  03 

1096  2  02  1146  XFR  67  1196  +  33 

1097  XFR  67  1147  X  36  1197  XFR  67 

1098  X  36  1148  3  03  1198  3  03 

1099  3  03  1149  XFR  67  1199  XFR  67 
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1200  X  36  1250    7  07  1300  DIV  35 

1201  4  04  1251    +  33  1301  1  01 

1202  XFR  67  1252  XFR  67  1302  0  00 

1203  X  36  1253    2  02  1303  XTO  23 

1204  6  06  1254  XFR  67  1304  2  02 

1205  +  33  1255    X  36  1305  XFR  67 

1206  XFR  67  1256    4  04  1306  1  01 

1207  3  03  1257  XFR  67  1307  3  03 

1208  XSQ  12  1258    X  36  1308  XFR  67 

1209  XFR  67  1259    8  10  1309  DIV  35 

1210  X  36  1260    +  33  1310  1  01 

1211  7  07  1261  XFR  67  1311  0  00 

1212  -  34  1262    3  03  1312  XTO  23 

1213  XFR  67  1263  XSQ  12  1313  3  03 

1214  2  02  1264  XFR  67  1314    .  21 

1215  XFR  67  1265    X  36  1315  5  05 

1216  X  36  1266    8  10  1316  XFR  67 

1217  5  05  1267    -  34  1317  DIV  35 

1218  XFR  67  1268  XFR  67  1318  1  01 

1219  X  36  1269    2  02  1319  CHS  32 

1220  6  06  1270  XFR  67  1320  UP  27 

1221  -  34  1271    X  36  1321  XFR  67 

1222  XFR  67  1272    5  05  1322  6  06 

1223  1  01  1273  XFR  67  1323  5  05 

1224  XFR  67  1274    X  36  1324  XSQ  12 

1225  X  36  1275    7  07  1325  X  36 

1226  4  04  1276    -  34  1326  YTO  40 

1227  XFR  67  1277  XFR  67  1327  2  02 

1228  X  36  1278    4  04  1328  5  05 

1229  8  10  1279  XSQ  12  1329  DN  25 

1230  /"  34  1280  XFR  67  1330  -  76 

1231  YTO  40  1281    X  36  1331  XTO  23 

1232  1  01  1282    6  06  1332  1  01 

1233  2  02  1283    -  34  1333  4  04 

1234  XFR  67  1284  YTO  40  1334  CNT  47 

1235  3  03  1285    1  01  1335  1  01 

1236  XFR  67  1286    3  03  1336  0  00 

1237  X  36  1287  XFR  67  1337    .  21 

1238  5  05  1288    1  01  1338  2  02 

1239  XFR  67  1289    1  01  1339  UP  27 

1240  X  36  1290  XFR  67  1340  5  05 

1241  6  06  1291  DIV  35  1341  FMT  42 

1242  UP  27  1292    1  01  1342  1  01 

1243  XFR  67  1293    0  00  1343  UP  27 

1244  3  03  1294  XTO  23  1344  FMT  42 

1245  XFR  67  1295    1  01  1345  1  01 

1246  X  36  1296  XFR  67  1346  FMT  42 

1247  4  04  1297    1  01  1347  D  63 

1248  XFR  67  1298    2  02  1348  A  62 

1249  X  36  1299  XFR  67  1349  XTO  23 
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1350  E  60 

1351  CNT  47 

1352  FMT  42 

1353  XFR  67 

1354  5  05 

1355  1  01 

1356  FMT  42 

1357  1  01 

1358  PNT  45 

1359  FMT  42 

1360  1  01 

1361  FMT  42 

1362  DIV  35 

1363  FMT  42 

1364  XFR  67 

1365  5  05 

1366  2  02 

1367  FMT  42 

1368  1  01 

1369  PNT  45 

1370  FMT  42 

1371  1  01 

1372  FMT  42 

1373  DIV  35 

1374  FMT  42 

1375  XFR  67 

1376  5  05 

1377  3  03 

1378  FMT  42 

1379  1  01 

1380  PNT  45 

1381  FMT  42 

1382  1  01 

1383  FMT  42 

1384  CNT  47 

1385  CNT  47 

1386  CNT  47 

1387  XTO  23 

1388  I  65 

1389  M  70 

1390  E  60 

1391  CNT  47 

1392  FMT  42 

1393  XFR  67 

1394  2  02 

1395  8  10 

1396  FMT  42 

1397  1  01 

1398  PNT  45 

1399  9  11 


1400 

. 

21 

1450 

1 

01 

1401 

9 

11 

1451 

PNT 

45 

1402 

UP 

27 

1452 

9 

11 

1403 

5 

05 

1453 

• 

21 
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6 
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27 
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UP 

27 
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5 

05 

1407 

FMT 

42 

1457 
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42 

1408 

1 

01 

1458 

1 

01 

1409 
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42 

1459 

UP 

27 

1410 

L 

72 

1460 

FMT 

42 

1411 

M 

70 

1461 

1 

01 

1412 

B 

66 

1462 

FMT 

42 

1413 

D 

63 

1463 

LBL 

51 

1414 

CNT 

47 

1464 

CNT 

47 

1415 

SFL 

54 

1465 

YTO 

40 

1416 

CNT 

47 

1466 

b 

14 

1417 

FMT 

42 

1467 

CNT 

47 

1418 

FMT 

42 

1468 

SFL 

54 

1419 

1 

01 

1469 

CNT 

47 

1420 

8 

10 

1470 

FMT 

42 

1421 

4 

04 

1471 

FMT 

42 

1422 

XFR 

67 

1472 

1 

01 

1423 

4 

04 

1473 

8 

10 

1424 

1 

01 

1474 

3 

03 

1425 

FMT 

42 

1475 

XFR 

67 

1426 

1 

01 

1476 

2 

02 

1427 

PNT 

45 

1477 

5 

05 

1428 

FMT 

42 

1478 

FMT 

42 

1429 

1 

01 

1479 

1 

01 

1430 

8 

10 

1480 

PNT 

45 

1431 

0 

00 

1481 

FMT 

42 

1432 

FMT 

42 

1482 

1 

01 

1433 

1 

01 

1483 

FMT 

42 

1434 

FMT 

42 

1484 

CNT 

47 

1435 

CNT 

47 

1485 

CNT 

47 

1436 

CNT 

47 

1486 

CNT 

47 

1437 

CNT 

47 

1487 

LBL 

51 

1438 

a 

13 

1488 

CNT 

47 

1439 

N 

73 

1489 

SFL 

54 

1440 

G 

15 

1490 

CNT 

47 

1441 

E 

60 

1491 

FMT 

42 

1442 

CNT 

47 

1492 

XFR 

67 

1443 

SFL 

54 

1493 

1 

01 

1444 

CNT 

47 

1494 

4 

04 

1445 

FMT 

42 

1495 

FMT 

42 

1446 

XFR 

67 

1496 

1 

01 

1447 

3 

03 

1497 
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45 
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1 

01 
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XFR 

67 

1449 

FMT 

42 

1499 

3 

03 

67 


1500  0  00  1550  FMT  42  1600  CNT  47 

1501  UP  27  1551    1  01  1601  INT  64 

1502  2  02  1552  FMT  42  1602  CNT  47 

1503  DIV  35  1553    C  61  1603  SFL  54 

1504  tt  56  1554    N  73  1604  CNT  47 

1505  DIV  35  1555  YTO  40  1605  FMT  42 

1506  7  07  1556    b  14  1606  FMT  42 

1507  UP  27  1557  CNT  47  1607  1  01 

1508  6  06  1558  SFL  54  1608  8  10 

1509  DIV  35  1559  CNT  47  1609  3  03 

1510  DN  25  1560  FMT  42  1610  XFR  67 

1511  XEY  30  1561  XFR  67  1611  1  01 

1512  H  74  1562    2  02  1612  CHS  32 

1513  XTO  23  1563    5  05  1613  1/X  17 

1514  X  36  1564  GTO  44  1614  XFR  67 

1515  2  02  1565  S/R  77  1615  X  36 

1516  5  05  1566  LBL  51  1616  2  02 

1517  XFR  67  1567    G  15  1617  UP  27 

1518  3  03  1568  FMT  42  1618  2  02 

1519  1  01  1569    1  01  16.19  DIV  35 

1520  UP  27  1570  FMT  42  1620  XFR  67 

1521  1  01  1571  CNT  47  1621  2  02 

1522  1  01  1572  CNT  47  1622  3  03 

1523  UP  27  1573  CNT  47  1623  -  34 

1524  6  06  1574    C  61  1624  XFR  67 

1525  DIV  35  1575    N  73  1625  6  06 

1526  DN  25  1576  CNT  47  1626  5  05 

1527  XEY  30  1577  SFL  54  1627  X  36 

1528  H  74  1578  CNT  47  1628  XFR  67 

1529  XTO  23  1579  FMT  42  1629  7  07 

1530  DIV  35  1580  XFR  67  1630  9  11 

1531  2  02  1581    2  02  1631  I  65 

1532  5  05  1582    5  05  1632  +  33 

1533  .  21  1583  /"  76  1633  XFR  67 

1534  4  04  1584  GTO  44  1634  1  01 

1535  9  11  1585  S/R  77  1635  4  04 

1536  6  06  1586  LBL  51  1636  XSQ  12 

1537  XTO  23  1587    G  15  1637  UP  27 

1538  DIV  35  1588    9  11  1638  2  02 

1539  2  02  1589    UP  27  1639  DIV  35 

1540  5  05  1590    5  05  1640  DN  25 

1541  CNT  47  1591  FMT  42  1641  +  33 

1542  9  11  1592    1  01  1642  DN  25 

1543  .  21  1593    UP  27  1643  J  75 

1544  3  03  1594  FMT  42  1644  XTO  23 

1545  UP  27  1595    1  01  1645  5  05 

1546  5  05  1596  FMT  42  1646  7  07 

1547  FMT  42  1597    A  62  1647  FMT  42 

1548  1  01  1598  INT  64  1648  1  01 

1549  UP  27  1599    E  60  1649  PNT  45 
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E 
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1759 

XEY 

30 

1660 
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70 
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CNT 
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SFL 

54 

1713 

UP 
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LBL 

51 
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5 

05 

1764 

8 

10 

1665 

FMT 

42 
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